TO APPEAR IN ASTROPHYSICAL JOURNAL 
Preprint typeset using I^TgX style emulateapj v. 10/09/06 



AEGIS: THE NATURE OF THE HOST GALAXIES OF LOW-IONIZATION OUTFLOWS AT Z < 0.6 



Taro Sato 12 , Crystal L. Martin 2 3 , 



Kai G. Noeske 4 , David C. Koo 5 , Jennifer M. Lotz 6 



(Received April 27, 2008; Accepted February 9, 2009) 

Draft version April 20, 2009 



o 
o 

u 

Oh: 

< 

o 

6 
& 



ABSTRACT 

We report on a signal-to-noise (S/N) limited search for low-ionization gas outflows in the spectra of the 
0.11 < z < 0.54 objects in the Extended Groth Strip (EGS) portion of the Deep Extragalactic Evolutionary 
Probe 2 (DEEP2) survey. Doppler shifts from the host galaxy redshifts are systematically searched for in the 
Nal A 5890,96 doublet (Nal D). Although the spectral resolution and S/N limit us to study the interstellar gas 
kinematics from fitting a single doublet component to each observed Na I D profile, the typical outflow often 
seen in local luminous-infrared galaxies (LIRGs) should be detected at > 6cr in absorption equivalent width 
down to the survey limiting S/N (~ 5 pixel" 1 ) in the continuum around Nal D. The detection rate of LIRG- 
like outflow clearly shows an increasing trend with star-forming activity and infrared luminosity. However, by 
virtue of not selecting our sample on star formation, we also find a majority of outflows in galaxies on the red 
sequence in the rest-frame (U-B, Mb) color-magnitude diagram. Most of these red-sequence galaxies hosting 
outflows are of early-type morphology and show the sign of recent star formation in their UV-optical colors; 
some show enhanced Balmer H/3 absorption lines indicative of poststarburst as well as high dust extinction. 
These findings demonstrate that outflows outlive starbursts and suggest that galactic-scale outflows play a role 
in quenching star formation in the host galaxies on their way to the red sequence. The fate of relic winds, 
as well as the observational constraints on gaseous feedback models, may be studied in galaxies during their 
poststarburst phase. We also note the presence of inflow candidates in red, early-type galaxies, some with signs 
of active galactic nuclei/LINERs but little evidence for star formation. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — galaxies: stellar content — 
ISM: jets and outflows 
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1. INTRODUCTION 

A number of large-scale galaxy surveys of unprecedented 
volume and depth are steadily mapping out the luminous con- 
stituents of the universe out to higher redshifts. The lambda 
cold dark matter (ACDM) cosmological paradigm describes 
the gravity-driven assembly history of the universe quite suc- 
cessfully on the largest physical scales. Yet our understand- 
ing of the universe remains largely incomplete on the galaxy 
scales where the astrophysical processes involving stellar evo- 
lution and feedback of interstellar and intergalactic matter 
are essential — baryons can cycle through these in an intri- 
cate, multiphase manner. Gaseous feedback processes have 
been gaining serious attention, since theoretical predictions 
of the growth of luminous structure are highly sensitive to the 
"prescriptions" for these small-scale, subgrid physics (e.£ 
ICroton et all 120061; iBower et alj 120061; iHopkins et all 12006 
yet the observational constraints are still relatively scarce. 

Absorption-line probes of the gas entrained in galactic 
"supe rwinds" (e.g.. lChevalier & Clegdll985l : lHeckman et alj 
1990) have been effective in constraining how much matter 
could be carried by outflows through the lines in the rest- 
frame o ptical (e.g.. | Heckman et all l2000t iRupke et al.1 12002), 
I2005bt iMartinl |2005j) and ultraviolet 
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1200 It ISchwartz et alJl2006l) . which may pollute intergalac- 
tic media and eventually deplete host galaxies of fuel for 
further star formation. The low-ionization absorption line 
studies suggest that a considerable amount of neutral gas 
may be carried away by o utflows, ~ 1 4 — 1 7 M Q in dwarfs 
(ISchwartz & Martinl l2004) and ~ lOMO 1 " in ultralumi- 
nous in frared galaxies (ULIRGs; Martin 2005; Rup ke etaT] 
2005b). The column density of outflowing cool matter, how- 
ever, can only be estimated after making a series of crude 
assumptions on its geometry, ionization state, dust-depletion 
factor, and metallicity. Although the reliability of such mass 
outflow estimates could therefore be questioned, an unam- 
biguous detection of outflow is relatively straightforward in 
finding a Doppler-shifted absorption component in a line 
complex of interest. So far, a vast majority of existing studies 
of outflows have focused on galaxies selected a priori to have 
high star formation rates (SFRs), fro m dwarf starburst galax- 
ies (e.g., ISchwa rtz & Martin 200 4]) and luminous infrared 
galaxies re.g.. lMartinl2005 . 2006; Rupk e et alJ2005allbl) in the 
local universe to Lyman-b reak galaxies at hig h redshifts (e.g., 
Shapl ev et alj|2003l) : see lVeilleux etail d2005l) for a recent re- 
view. Although winds have been unambiguously detected in 
these systems, much work remains to be done in characteriz- 
ing outflows in terms of their host galaxy properties in a large, 
relatively unbiased sample. 

The Na I A 5890, 96 doublet (Na I D) absorption can be stel- 
lar or interstellar in origin yet is a useful line for a census of 
cool winds in the interstellar medium (ISM). Since it is an ab- 
sorption line measured against the background light of host 
galaxies, its Doppler shift relative to the systemic redshift can 
be measured reasonably well. The Na I D absorption line di- 
rectly traces cool (T ~ 100 K) gas which may directly fuel 
star formation. In (U)LIRGs, ~ 10% of the dynamical mass 
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can b e entrained in cool outflows (Rupke et al. 2005b; Martin 
2006), suggesting that the bulk of outflowing mass resides at 
the cool phase. The Na I D doublet is a resonance line and 
forms among the most prominent absorption lines detectable 
out to z ~ 0.5 in optical spectra. 

There is a pressing need to extend the study of outflows to 
the systems with lower SFRs as well as at later stages of star 
formation. For one thing, the observed spatial extent of out- 
flows in local ULIRGs and a simple dynamical argument sug- 
gest that outflowing gas clouds may outlive starbursts or ac- 
tive g alactic nucle i (AGNs) which may dr ive outflows (Martin 
120061) . Recently, iTremonti et aD (120071) found that z ~ 0.5 
massive poststarburst galaxies, observed at up to 1 .5 Gyr after 
intense episodes of star formation, host outflows almost as fast 
as or even faster than those observed in local ULIRGs. Since 
outflows may carry away from galaxies a significant fraction 
of gas mass that would otherwise be available for further star 
formation, knowing the "fate" of outflowing matter is clearly 
of interest. 

In this paper, we present the result from a systematic 
search for Na I D outflows in a flux-limited sample of galax- 
ies, drawn from the Extended Groth Strip (EGS) portion of 
the Deep Extrag alactic Evolutionary Probe 2 survey (DEEP2; 
iDavis et alj|2003l) . Since the EGS field has been extensively 
observed by a wide array of multiwavelength missions as part 
of the All-wavele ngth Extended Groth Strip Survey (AEGIS; 
IDavis etalJl2007l) . we are able to characterize our spectro- 
scopically selected sample in view of various physical quan- 
tities. Although Na I D is a resonance line and among the 
most prominent of absorption lines in the visible, rigorous ab- 
sorption analysis needs high signal-to-noise (S/N) continuum, 
which is a requirement not quite satisfied by the vast majority 
of 1 hr integration spectra from the DEEP2 survey. Co-adding 
a number of low-S/N spectra from a set of subsamples to im- 
prove the effective S/N is a sensible approach. Yet even such 
stacking analyses are limited a priori by numerous ways in 
which subsamples can be constructed; some prior knowledge 
must be gained about the subsampling schemes in which the 
desired information can best be elucidated. In order to initi- 
ate an effort to carry out an unbiased census of outflows in 
modern spectroscopic galaxy surveys, as well as to motivate 
ensuing stacking analysis with proper subsampling schemes, 
it is still beneficial to take an approach to seek evidence of 
outflows in individual spectra. 

An ambitious goal would be to estimate the quantities that 
are useful for improving the prescriptions of cosmological 
semianalytical simulations, such as outflow detection rate and 
mass-loading factor in a robust volume-limited sample. In 
this paper, however, our focus is on simply characterizing 
the property of galaxies that host LIRG-like outflows, using 
a rich set of multiwavelength observations from the AEGIS 
survey. We first describe the data, selection, and analysis 
method for the Na I D sample in §|2] The host galaxies of out- 
flows, as defined from the Na I D kinematics, are then studied 
in view of their UV, optical, and infrared properties in §[3] 
We then note a few caveats, frame our findings in the larger 
context of galaxy evolution (§0), an d summarize the paper 
(§0). Throughout this paper, we adopt the standard ACDM 
cosmology, (fi m , A ) = (0.3, 0.7), with H Q = 70 km s" 1 Mpc" 1 . 

2. DATA AND ANALYSIS 

2.1. DEEP2 Spectra 



The EGS, covering w 0.5 deg 2 , is one of the four fields ob- 
served in th e DEEP2 survey (see lCoil et alJl2004tlDavis et alJ 
I2003L 12004 for the descriptions of the survey), in which the 
spectroscopic targets are preselected by the Canada-France- 
Hawa ii Telescope (CFHT) BRI photometry (ICuillandre et alJ 
1200 ll) . While the targets in the other three fields are pres- 
elected to the apparent magnitude limit of /?ab < 24.1 and 
by color cuts to the galaxies likely to be at z > 0.7, such a 
color cut was used in the EGS to slightly down-weight galax- 
ies at z < 0.7, resulting by design in roughly equal num- 
bers of galaxies above and below z = 0.7. The spectra for 
~ 13,570 EGS objects are obtained with the DEIMOS spec- 
trograph on the Keck II telescope ( Fa ber et all 12003 ) . The 
1200 lines mm" 1 grating centered at 7800 A and the OG550 
order-blocking filter with 1" slit widths are used, which leads 
to the spectral coverage roughly of 6500-9100 A. The spec- 
tral resolution in FWHM is ~ 68 km s" 1 . The data are pro- 
cessed by an automated pipeline to produce (unfluxed) one- 
dimensional spectra (Newman et al., in preparation). The 
DEEP2 DEIMOS pipeline 7 produces an inverse-variance vec- 
tor for each spectrum. In this paper, we use the spectra ex- 
tracted vi a a var i ant of the optimal extraction algorithm pre- 
sented in lHornd d!986l) . Throughout this paper, we only use 
spectra which are tagged DEEP2 ZQUALITY flag of four. 

2.2. Systemic Redshifts 

An outflow velocity is measured relative to a systemic red- 
shift of a host galaxy. A major contribution to the uncer- 
tainty of outflow velocity therefore comes from the uncer- 
tainty in the systemic redshift measurement. Each DEEP2 
spectrum is assigned a spectroscopic redshift from \ 2 min- 
imization with a few template spectra and has the accuracy 
well quantified from repeat observations (30 km s" 1 in RMS; 
IWillmer et al.ll2006h . Each redshift measurement has been vi- 
sually inspected and assigned a redshift quality flag. 8 Never- 
theless, we carry out independent reds hift measurements us - 
ing IRAF software package XCSAO (iKurtz & Minklll998l) . 
The primary motivation for an independent set of redshift 
measurements is to mask the Na I D complex. The Na I D 
absorption line is among the most prominent features in the 
visible part of the galaxy spectrum. In the case where a sys- 
temic redshift should only reflect the centroid of the stellar 
motions, Na I D should be excluded from cross correlation 
with stellar spectral templates, in order to avoid the interstel- 
lar components, which may be redshifted or blueshifted from 
the systemic redshift of the galaxy, to affect the result. 

We adapt the cross-correlation templates (ID: 24-30) from 
the fifth data release of the Sloa n Digital Sky Survey (SDSS; 
lAdelman-McCarfhv et alj|2006l) . The template which yields 
the smallest uncertainty is generally picked and designated as 
the systemic redshift of the galaxy. An uncertainty in each 
redshift measurement is derived from the r statistics, which 
roughly corresponds to the S/N of cross-correl ation peak from 
which the best redshift estimate is computed (Tonrv & Davis 
119791) . The agreement between XCSAO and DEEP2 red- 
shifts is generally good; for most purposes, the systematic 
redshift disparity of cz(XCS AO) - cz(DEEP2) » 10 km s" 1 is 
much less than the instrumental resolution and is insignificant. 
Where precise systemic redshifts are critical, we use XCSAO 
measurements, since the uncertainty is empirically calibrated 
and well understood among our sample. The redshift mea- 

7 http://astro.berkeley.edu/~cooper/deep/spec2d/ 

8 See http://deep.berkeley.edu/ for the DEEP2 survey detail. 
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surement stored in each DEEP2 spectrum almost certainly 
underestimates the uncertainty in case of gross template mis- 
match, since only three components are used to generate a 
template for cross correlation, and \ 2 statistics was computed 
after continua are removed from them. Each redshift measure- 
ment for which |cz(XCSAO)-cz(DEEP2)| > 2o\,(XCSAO) is 
visually inspected, and the spectrum is either reassigned a red- 
shift from another template or removed from the sample. Only 
one object dropped out of the sample because of the failure in 
the cross-correlation redshift measurement. 

2.2.1. Spectral Line Indices 

Using the rest-frame band definitions in TableQ] we mea- 
sure the spectral indices (denoted Wo, with the subscript "0" 
meaning rest frame) of some lines as estimates of their equiv- 
alent widths. Each spectral line index is computed via the 
"flux-summing" method. First, a straight line "pseudocontin- 
uum" is fitted to the variance-weighted pixel flux values from 
the blue and red straddling continua as defined in Table Q] via 
the standard Levenberg-Marquardt algorithm. The covariance 
matrix is used to estimate the uncertainty in the pseudocon- 
tinuum. Then at each pixel i with a pixel width AA,- in wave- 
length falling within the spectral bandpass defined as "line" 
in TableQ] the observed flux f and the pseudocontinuum flux 
f c is used to compute the flux excess or depletion, such that 



TABLE 1 
Spectral Line Index Definitions 



W 



-E 



fo-fc 
fc 



AX, 



yields the equivalent width index in the observed frame. The 
rest-frame value is computed simply from Wo = W/(l+z) 
where z is the redshift of the galaxy. The uncertainty in the 
index is estimated by formal propagation of uncertainties in 
f and f c using the above relation. By construction, Wo > A 
(< A) corresponds to a line flux seen in absorption (emis- 
sion), although the physical interpretation is slightly compli- 
cated by the fact that both emission and absorption may exist 
in a line feature, sometimes called emission filling. Although 
not perfect, a line index gives a good estimate of true equiva- 
lent width, when emission filling is not significant. 

Previous studies have used the equivalent width of the 
Mg I bA5167,73,85 triplet to estimate the stellar contribu- 
tion to the Nai D absorption line (|Heckman et al.l 120001; 
Rupke etalJ 120021: [Schwartz & Martini 12004 iRupke et al. 
2005ai: Martin 2005). The correlation between NalD and 



Mg I b equivalent widths in stellar spectra is expected based 
on the similar mechanisms from which Na and Mg are pro- 
duced in stars and their roughly similar ionization poten- 
tials (5.14 eV a nd 7.6 5 eV for Na and Mg, re spectively). 
HeckmanetaD d2000l) . iMartinl (120051) . and IRupke et all 
(I2005ah all showed high equivalent width ratios of Na I D to 
Mg I b to be a good indicator of the presence of winds from 
their samples of infrared-luminous galaxies. This assumption 
is reasonable provided that a presence of a large column den- 
sity from the interstellar Na I is required for a secure detection 
of outflow. In Fig.Q] however, we see evidence that a pop- 
ulation of galaxies with outflows would be missed by a high 
W)(Na I D)/Wo(Mg I b) selection scheme. The figure is meant 
to facilitate a comparison to the existing studies, in which the 
samples are selected by a high level of star formation (i.e., a 
high infrared luminosity). As will be discussed later, outflows 
appear to outlive starbursts, and the relic winds may present 
significant columns to be detected in poststarburst or post- 
star-forming galaxies. 
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2.3. Sample Selection 

We systematically search DEEP2 spectra in the EGS field 
for the coverage of Na I D absorption line and continuum 
around it (TableQ}. An object is removed from our sample 
if the redshifted Na I D spectral range is at least partially out- 
side the edge or falls on the gap between the blue and red CCD 
chips. This process reduces the sample to 2248 objects. The 
spectral baseline of DEEP2 observations restricts our sample 
to 0.11 < z < 0.54. 9 For each DEEP2 spectrum, the contin- 
uum S/N per pixel [(S/N) p i x ] in the region around Nai D is 
defined to be the median (S/N) p j x computed from the pixel 
flux values and their inverse variances registered within the 
blue and red continua (TableQ}. The main selection cut is 
made at (S/N) p ; x > 5 around Na I D, reducing the sample to 
493 objects. We also remove the objects whose Na I D fea- 
ture is severely compromised by sky emission or atmospheric 
absorption lines after visual inspection, reducing the sample 
to 431 objects. The latter cut tends to remove objects at spe- 
cific redshifts where the redshifted Na I D overlaps with tel- 
luric features. 

The particular choice of continuum S/N cut is a compro- 
mise between the inclusion of more objects for better statistics 
and the reliability of Na I D velocity measurements. Given the 
limited spectroscopic S/N, our desire to probe fainter objects 
is motivated by the well-known , downsizing nature of star for- 
mation (e.g. JCowie et al.ll 19961) and the apparent correlation 
between the presence of galactic-scale wind and the stre ngth 
of star formation (e.g.. lMarfirJl2005t IRupke et al.ll2005bl) . At 
low z, star formation is expected in optically fainter, low 
surface brightness galaxies. At high z, increasingly brighter 
galaxies are host to star formation yet become faint in their 
apparent brightness in the visible. We also expect that the 
high-S/N spectra are obtained from passive, early-type galax- 
ies, which are generally of high surface brightness. All these 
effects conspire to make the galaxy population of interest to 
be somewhat elusive in the optical selection used here. Fur- 
thermore, the strength of Na I D absorption, both stellar and 
interstellar, varies widely from galaxy to galaxy, and our abil- 
ity to detect an outflow depends on the strength of continuum 
as well as the absorption feature. It is therefore important to 
get some idea as to what kind of Na I D outflow to which our 
measurements are sensitive in this study. 

In Fig. |2] the distribution of spectral index measurements as 
a function of continuum (S/N) p i x around Na I D is shown and 
demonstrates that we cannot obtain reliable Na I D velocity 
measurements in most spectra (i.e., constitutes the low-S/N 
sample; see §Q3]for detail) at (S/N) pix < 6.5. The detec- 
tion limit as a function of (S/N) p j x for a fiducial LIRG wind 
indicate that we are reaching the ~ 6a detection limit for a 
typical LIRG-type wind in our sample at that continuum S/N 
level. The distribution of low-S/N measurements suggests 
that we lose the ability to detect and measure the kinemat- 

9 Over the redshift range, the physical scale corresponding to 1" slit varies 
from 2 kpc to 9.5 kpc. 
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FIG. 1 . — Mg I b-Nal D spectral index plane for the objects with the measurements of Na I D velocity (left) and the infrared luminosity (right). See § 12.2.11 
for the definition of spectral indices. In the left panel , the p oints are denoted by Nal D kinematics into outflow (blue filled circle), systemic (gray cross), inflow 
(red open cir cle), an d low S/N (light gray dot); see § 12.4. 1 I for the detail on Nal D kinematics. Outflows/inflow data points are color-mapped with the blueshift 
probability (§ 12.4. It . where bluer (redder) marks indicate stronger blueshift (redshift) of Na I D absorption line; see the inset of Fig.[5]for the exact color-mapping 
scheme with the blueshift probability. The error bars shown for the spectral indices are the medians for high- and low-S/N samples. In the rig ht panel, the 
symbol s indicate LIRGs (red open circle), non-LIRGs (blue crosses), and objects below flux limit of the Spitzer/MIPS 24 pm survey (gray dot). IRupke et all 
1 2005a) found most of their infrared- selected objects lying below the gray dotted line, W(NaI D) = 3W(Mg I b), to host outflows, whereas most objects above 
did not. The long gray dashed line indicates the fiducial stellar loc i, W(Na I D) = 0.4iy (Mg I b), where the slope was chosen to roughly match the evolutionary 
track of the solar-metallicity, single stellar population model from Delgado et al. 1 2005, solid green line). For each object the differences between the observed 
and fiducial stellar Nal D indices presumably indicate the interstellar contribution (see Fig.[3J. (A color version of this figure is available in the online journal.) 



ics of LIRG-like outflows at (S/N) p ; x w 5. Thus the selec- 
tion cut at (S/N)pix = 5 seems justified, in terms of detecting 
LIRG-like winds at high (> 5a) confidence. Fig.[2]shows that 
the success rate of Na I D velocity measurement is a strong 
function of continuum (S/N) p ; x . Overall, there are 205 ob- 
jects with successful Na I D velocity measurements, and 226 
objects without. We will describe what constit utes a " success- 
ful" (i.e., high S/N) velocity measurement in § 12.4.11 



2.4. Modeling Na I D Absorption Lines 

In principle, countless possible configurations of the ge- 
ometry of individual absorbers along a sightline give rise to 
an unlimited variety of observed Na I D absorption line pro- 
files; from numerical simulations, several absorbers entrained 
in a starburst wind are expected to lie along a single sight 
line (A. Fujita et al. 2008, in preparation). In practice, how- 
ever, the spectral resolution and moderate S/N limit our abil- 
ity to study more than one component of Na I D doublet in 
the DEEP2 spectra; multiple absorption components would 
be seen blended even at a sufficient S/N. Due to the limita- 
tion, the physical quantities derived from a Na I D line profile 
in general may not be uniquely determined. Nonetheless, the 
modeling of the absorption line should be physically moti- 
vated, an d in this resp e ct we closel y follow the method pre- 
sented by Rup ke et"all d2002l l2005al) ; readers are highly en- 
couraged to find the detail of their absorption line analysis 
method in those papers. An absorption line is modeled with 
the wavelength (A c ) and the optical depth (to) at the line cen- 
ter, Doppler width {bo), and the covering fraction (C/) in a 



self-consistent manner. Our confidence intervals on Na I D 
velocity measurements, however, are obtained via directly 
carrying out the Markov chain Monte Carlo (M CMC) analysis 
on the observed spectra wit h the | Rupke et all- typ e line-profile 
modeling, rather than what IRupke et al.l (l2005al) outlines. In 
the Appendix, we give a summary of the technique, the anal- 
ysis method, and the pipeline software developed for the task. 

MCMC sampling generates a probability density for each 
model parameter which allows us to visually inspect the qual- 
ity of our measurements. Except for a couple dozen high S/N 
spectra, the optical depths cannot be constrained at all, i.e., 
the probability density for the central optical depth tq usu- 
ally ends up being distributed roughly uniformly over the al- 
lowed range (0 < to < 999) with a slight enhancement to- 
ward the lowest optical depths. In turn, the highly saturated 
profile tends to let the covering fraction C/ be distributed 
near the level of minimum intensity of a Na I D profile. This 
degeneracy between to and Cf at a low S/N regime is a 
well-unders tood property of the model profile employed by 
Rupk e etaLl and in this paper; extracting an optical depth in 
general requires a line shape to be very well sampled. The 
distributions of Na I D central wavelengths A c and Doppler 
widths bo, on the other hand, are relatively well behaving 
even at lower S/Ns, where their probability densities roughly 
become Gaussian. 

One especially important caveat of our Na I D velocity mea- 
surement is in order. Our interest is in studying the interstel- 
lar gas kinematics, so ideally stellar contributions to Na I D 
should be removed via such a method as fitting template 
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Continuum (S/N) p ; x around Na I D 

FIG. 2. — S/N (top) and spectral line index (bottom) as a function of con- 
tinuum (S/N) p j x around Na I D. The points are denoted by Na I D kinematics: 
outflow (blue filled circle), system ic (gra y cross), inflow (red open circle), 
and low S/N (light gray dot); see 8 12.4. ll for the detail on Nal D kinematics. 
O utflow s/inflow data points are color-mapped with the blueshift probability 
(§ 12.4. U . where bluer (redder) marks indicate stronger blueshift (redshift) of 
Nal D absorption line. The green lines are the loci of Nal D indices mea- 
sured in the simulated DEEP2 spectra assuming fiducial Na I D absorption 
profiles expected from the fixed optical depth and covering fraction (t = 1.1 
and Cf = 0.4) at three different Doppler widths [bo = 50 (dotted green line), 
150 (dashed green line), and 350 km s _1 (dash-dotted green line)], which 
are varied to simulate the effect of broadening in the absorption profile (not 
convolved with the instrumental resolution of rs 42 km s -1 here); see the Ap- 
pendix for the definitions of the variables. The model with bo = 150 km s~* 
roughly reflects the mean Na I D property of LIRGs reported by Ru pke et"aT] 
1 2005b); i.e., the long dashed line in the upper panel indicates the significance 
of detection for a typical LIRG wind. For each model, a number of realiza- 
tions are generated at a given S/N (i.e., a spectrum is degraded by Gaussian 
noise assuming that S/N) to obtain the uncertainty in the spectral index. (A 
color version of this figure is available in the online journal.) 

spectra generated fr om population synthesis models (e.g., 
iTremonti et al.l l2004). The DEEP2 spectra are not rigorously 
fluxed, however, and we are unable to carry out a similar 
procedure. This is a substantial limitation in the analysis of 
galaxy spectra of intermediate- and old-age stellar popula- 
tions, since their stellar absorpti on at Na I D becomes strong 
(Fig. [j] the stellar loci are from lDelgado etal.1 d2005l) ). The 
implication on our definition of outflow velocity will be dis- 
cussed in § 12.51 

We also do not take any special care of the nebular emis- 
sion line He I A 5876, found ~ 15 A blueward of Nal D in 
some spectra. The He I emission line can contaminate the 
high- velocity tail of strong Na I D outflows in (U)LIRGs (e.g., 
Rupk e et alTl2005at iMartinl 120051) . Upon visual inspection, 
however, few galaxies in our sample with a good Na I D ve- 
locity measurement are found to have their Na I D contami- 
nated by the presence of the He I emission line. 

2.4. 1 . Na I D Velocity and Blueshift Probability 

Each Na I D central wavelength A c is converted to the 
Na I D velocity offset 



v(Na I D) = c 



A c — A s < 
Asys 



(1) 



where c is the speed of light and A sys = 5895.9243(1 + z) A 
is the line center of Na I D shifted to the observed frame us- 
ing the cross-correlation redshift z. (The redder line of the 
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FIG. 3. — Distribution of Nal D velocity confidence intervals as a function 
of continuum (S/N) p ; x around Na I D. Here Av^% is the difference of upper 
and lower 68% confidence limits in Na I D velocity. The symbols are as in 
Fig. [2 Inset: the Nal D velocity measurement with the error bar indicating 
the 68% confidence interval as a function of blueshift probability pbiue- Out- 
flows and inflows from the high-S/N sample are color-coded by blue and red, 
respectively, in gradient according to Pbiuc (sectNalD Velocity and Blueshift 
Probabilit y); tho se at the systemic velocity are in gray. The low-S/N velocity 
sample (§ 12.4. U are in light gray. High- and low-S/N velocity measurements 
are effectively separated in the ways that they are distributed in these plots, 
indicating the efficacy of the visual inspection scheme; see the text for detail. 
(A color version of this figure is available in the online journal. ) 

doublet, Na I A 5896, is used as the reference line through- 
out.) Using the above equation, the probability distribution 
of v(Na I D) is directly obtained from the probability distri- 
bution of A e from the MCMC sampling for each spectrum. 
In order to take into account the systemic redshift uncertainty, 
the probability distribution is further convolved by a Gaussian 
kernel having a width corresponding to the 1 a uncertainty in 
the cross-correlation redshift (in the velocity space) for each 
object. From each probability distribution, the best estimate 
for Nal D velocity is taken from the median, and the 68% 
confidence interval is likewise obtained. 

For the purpose of measuring the Na I D velocity, the ab- 
sorption feature detected at a sufficient S/N for such a mea- 
surement generally reveals itself as a normally distributed 
probability density in v(Na I D) well bounded within the pa- 
rameter range, v(Na I D) = ±700 km s" 1 . 10 At a lower S/N, a 
probability distribution for v(Na I D) exhibits high- and low- 
velocity tails reaching the boundary values. Hence, after vi- 
sual inspection, we divide our v(Na I D) measurements into 
two classes based on the behavior of v(Na I D) probability 
distribution: "high-S/N" velocity sample (TV = 205) for which 
the distribution is well within ±700 km s" 1 and "low-S/N" 
velocity sample {N = 226) for which the distribution either 
extends to the boundary or is ill-behaving. Visual inspection 
also guards against the sampling results latching on to un- 
wanted noise features, which usually show up as an abnormal 
probability distribution function. The MCMC measurement 
pipeline also allows a fitted absorption-line profile to be in- 
spected for an interactively picked set of model parameters, so 
the integrity of the fitting result has also been visually checked 

10 While the Nal D outflow velocities outside this int erval have been re- 
ported in literature (e.g., Rupke et al. 2005b; Martin 2005), the visual inspec- 
tion of the velocity measurements indicates that the range of ±700 km s - ' 
is sufficient for our sample, which does not include ULIRGs (with signs of 
AGN activities) in which the most very high velocity outflows are observed. 
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FIG. 4. — Top: examples of continuum-normalized spectra around Nal D 
(black line) and the best-fit models for the blueshift probability Pblue — 1 
(blue), 0.5 (gray), and 0.01 (red). The vertical gray dashed line indicates 
the systemic velocity. The vertical dotted lines indicate the line center of 
Nal A 5896 from the best-fit model for each spectrum. The km s cor- 
responds to the systemic velocity of the red line of a blended doublet, and 
thus appears shifted toward red with respect to the centroid of a line profile 
even if the doublet itself is at systemic. Bottom: the p robabi lity distributions 
of Nal D velocity from the MCMC sampling; see 8 12.4. II for the detail on 
how the distributions are obtained. For the object showing a very strong out- 
flow (blue), the probability distribution is entirely found at < km s , so the 
blueshift probability />bi U c approaches unity. For a majority of galaxies, the 
distributions are roughly symmetric about km s _1 , for which Pblue is about 
0.5 (black). (A color version of this figure is available in the online journal.) 

at various points in the distributions of model parameters; see 
the Appendix for detail. Fig. [3] shows the velocity width of 
confidence intervals (Av^%; i.e., the difference between the 
upper and lower 68% confidence intervals) as a function of 
Nal D continuum (S/N) p ; x . Since Av'68% is a measure of the 
width of the probability distribution for v(Na I D), the figure 
shows that our strategy effectively distinguishes and separates 
out high- and low-S/N velocity measurements, which tend to 
appear at low and high Av68% regions, respectively. 

Since we have robust probability density functions for 
v(Na I D), it is desirable to incorporate these into our def- 
inition of outflow detection, rather than relying heavily on 
the best estimates for v(Na I D). To do this, we introduce 
the blueshift probability pblue for each velocity measurement, 
which is simply given by 

,o 

Pblue = / dvp(y) , 

J — oo 

where p(v) is the probability distribution function obtained 
from the MCMC sampling; recall that a blueshift yields a neg- 
ative velocity in Eq. (Q~|i, so the distribution function needs to 
be integrated out to negative infinity to obtain the probability 
that a Na I D is seen blueshifted in a given spectrum. A few 
examples of observed Na I D absorption spectra, along with 
the probability distributions of Na I D velocity, are shown in 
Fig- H] for the cases of pbiue — 1 (certainly an outflow), 0.5 
(likely at systemic), and 0.01 (almost certainly "inflow"). The 
value of pblue is a measure of how likely that a "real" Na I D 
velocity is blueshifted from the systemic velocity, given the 
result of MCMC sampling and the uncertainty in the systemic 
redshift of a host galaxy. 

2.5. Definition of Outflow 



To define what constitutes an outflow ("inflow") detec- 
tion for the current analysis, we use a Na I D velocity cut of 
v(Nal D) < -50 km s" 1 (> +50 km s" 1 ) and a blueshift prob- 
ability cut of pbiue > 0.75 (< 0.25). Readers are cautioned 
that, based on instrumental effects and data quality, precise 
definitions of outflows necessarily vary in the literature. The 
specific choice for the velocity cutoff is primarily motivated 
by the typical v(Nai D) uncertainty of ~ 50 km s" 1 (Pig. [3j 
and the visual inspection of fitting results in relation to the 
values of v(Nal D) and pbi ue - Although a blueshift probabil- 
ity, by construction, should be an indicator of the likelihood 
of detecting an in/outflow at a desired confidence level, the 
additional velocity cutoff is used to guard against both the 
random and systematic errors in the redshift and the velocity 
measurements. For example, when the width of probability 
distribution function is narrow, i.e., Av68% is small, pbi ue be- 
comes more sensitive to the particular value of the systemic 
velocity. Upon visual inspection, we found that this often 
happened when a Na I D doublet line shape was less blended, 
and therefore the velocity moment could be determined more 
precisely. In such a case, an outflow as defined solely by a 
v(Na I D) distribution function looks dubious, having a ten- 
dency to be thrown off by an error in a single measurement of 
systemic redshift. Ideally, our systemic redshifts should also 
be scrutinized under rigorous MCMC analysis, which we did 
not carry out. For now, the additional velocity cut adequately 
achieves th e same goal. The s imilar velocity cutoff was em- 
ployed by Rup ke et ail (f2005b) to define outflow. 

While employing a more stringent pbi ue cut would yield a 
sample of outflows detected with higher confidence, we then 
seriously have to compromise our sample size. In order to 
take advantage of the large amount of data available in the 
AEGIS survey, our approach is to push the limit of detections 
down to an acceptably low confidence level. Hence our goal 
in this study is not to find individual cases of outflows or in- 
flows securely but to give ourselves some statistical power to 
characterize the general properties of host galaxies. This is a 
compromise that we opt to consciously make. 

As emphasized earlier, one caveat of the Na I D absorption 
analysis in this paper is our inability t o rem ove the stellar 
contribution in the Nal D absorption (§ 12.41 : a stellar frac- 
tion can only be estimated indirectly from such an index as 
Mg I b whose strength is known to correlate well with that of 
Na I D. Since the presence of interstellar absorber(s) is a nec- 
essary condition for outflows seen in absorption, some trend 
is expected to exist between the outflow detection rate and 
the interstellar fraction of the total Na I column. Fig.|5]shows 
that such a trend does exist, in which the objects with low 
Na I D stellar absorption fractions are more likely to host out- 
flows (that we can detect). Without the explicit removal of 
stellar absorption components, however, the figure does not 
necessarily imply the paucity of outflows in galaxies whose 
Na I D is dominated by the stellar contribution. It is important 
to note again that our census is only sensitive to fairly strong 
outflows of the kind expected in LIRGs (Fig.E). The k ind 
of weak outflows observed by ISchwartz & Martini (|2004|) in 
dwarf starbursts is certainly below the detection level of the 
present survey. Furthermore, since only one Na I D doublet 
component is fitted, the measured Na I D velocities are likely 
the lower limits to the kinematic component with the high- 
est outflow velocity; each Nal D velocity is sensitive to the 
"moment" of multiple absorption components, at least one 
of which is stellar in origin and should exist a t a systemic 
redshift. In high-S/N spectroscopy of (U)LIRGs, Rup ke et all 
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uncertainties to be very high (greater than 0.7). Bottom: the fraction of ob- 
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objects. (A color version of this figure is available in the online journal.) 



(2005a|) and lMartir] d2005l) . for example, could fit more than 
one component of Na I D absorption in some spectra, pro- 
vided that line profiles are not overly "smooth" due to blend- 
ing and instrumental smearing. Fig. ^indicates that our veloc- 
ity measurements appear sensitive to outflows in Na I D when 
> 50% of the absorption equivalent width is interstellar in ori- 
gin. Given these limitations, attaching a physical meaning to 
a Na I D velocity would be misleading, since it is only the 
shift in the moment of the absorption profile, not the isolated 
absorption from the outflowing/inflowing gas, to which our 
measurements are really sensitive. We thus avoid emphasis 
on the exact Na I D velocity values. 

2.5.1 . Reality of "Inflows " 

The sensible cuts in the Na I D velocity and blueshift proba- 
bility pbiue give rise to a population of galaxies with "inflows." 
While it would not be surprising to see an inflow from a sight- 
line through an interacting system, "inflows" in our sample 
are seen mostly in luminous, massive galaxies in the red se- 
quence (§ 13. 2\ . The presence of inflows in these objects is 
striking. Given that these "inflow" detections are almost ex- 
clusively in early-type galaxies presumably with relatively lit- 
tle interstellar gas, we speculate whether factors other than 
interstellar absorption can cause potential systematics. 

First, the precision of cross-correlation redshifts is gen- 
erally worse with templates dominated by absorption lines, 
from which the redshifts of galaxies with "inflows" are ob- 
tained. It is also difficult to visually inspect redshift systemat- 
ics in absorption-dominated spectra due to the lack of nar- 
row, high-S/N features. As mentioned in § 12.21 a system- 
atic velocity difference of order 10 km s" 1 also exists be- 
tween cross-correlation and DEEP2 redshifts, although this 



is small compared to inflow velocities. We do, however, take 
the reds hift un certainty into account in our definition of in- 
flows (§ 12.4. U . Second, it might be possible that some unac- 
counted metal absorption features redward of Na I D may be 
causing systematics, redshifting the single-component Na I D 
profile fit. The exp erimentati on with synthesis spe ctra from 
iBruzual & Charlotfd2003h and lDelgado et all d2005l) indicates 
that the degree of systematics would not be as strong as ob- 
served, assuming that Na I D and nearby metal lines are prop- 
erly modeled. Third, as a significant fraction of galaxies with 
inflows appear to have lenticular morphology (Tremonti et al. 
2007, private communication), another possibility is that a 
distinct kinematic component of stellar motion might be de- 
tected in these systems, perhaps from a faint disk. 

While some of the above listed concerns are equally valid 
for outflows, in following sections we show compelling evi- 
dence that the detections of outflows are physically associated 
with star formation (or AGNs/LINERs). Such a clear physi- 
cal connection cannot be made for inflows at present, and a 
preliminary investigation to explore the nature of the "inflow" 
population is underway. In radio ellipticals, neutral hydrogen 
is often seen in inflow, plausibly feeding the nuclear activ- 
ity (e.g.. Ivan Gorkom et al.l IT989h : we do find some consis- 
tency with this scenario in our inflows. If these inflows are 
indeed associated with such feeding of AGNs, our method 
could provide an effective means of identifying the massive, 
early-type galaxies going through the "maintenance-mode" of 
AGN feedback. Nonetheless, we mostly defer the discussions 
of inflows to future papers. 

2.6. On Selection Effects 

Dependence on Na I D absorption strength. Our ability to 
measure a Na I D velocity depends primarily on the combi- 
nation of the absorption and the continuum strengths, which 
introduces selection biases. Our main selection cut on contin- 
uum S/N around Na I D feature favors luminous, high surface 
brightness obje cts, w hich tend to be early-type galaxies on the 
red sequence (§ 13.21 ). Star-forming galaxies in the blue cloud 
are much fainter at z < 0.5 from which our sample is drawn. 
In Fig.Q] we see that high Wb(Nal D) objects tend to be old 
galaxies with a high stellar fraction or young galaxies with 
a high interstellar fraction in their Na I D absorption. Our 
survey lacks sensitivity to absorption lines in less luminous 
star-forming galaxies, where much of star formation occurs at 
Z < 0.5. Therefore, that a significant fr actio n of our outflow 
detections is in red-sequence galaxies (§ 13.21 ) is partially a se- 
lection effect. Nonetheless, our sample nicely complements 
the existing studies of Na I D outflows which have focused on 
starburst galaxies. 

Redshift dependence. Our sample consists of objects dis- 
tributed over a fairly large redshift interval of 0. 1 1 < z < 0.54, 
and the analysis suffers from the problem typically encoun- 
tered by a flux-limited survey. First, the observation becomes 
increasingly insensitive to fainter populations at higher red- 
shift. Second, given a fixed opening angle, the survey volume 
tends to be smaller at lower z, and the sample variance (i.e., 
cosmic variance as galaxy enthusiasts like to say) becomes 
a serious issue; the survey field of view changes by a fac- 
tor of 10 in physical area over the redshift interval. Third, 
the lookback time difference of several Gyr means that the 
galaxies experience a significant evolution over that redshift 
interval, so the lower- and higher-z populations characterized 
by one "fixed" physical property may not be of the same kind 
when looked at in view of other properties. 
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FIG. 6. — Stellar mass of galaxy as a function of redshift for the sample in 
this study. Each galaxy is coded by the Nal D kinematics (Fig. [5): outflow 
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S/N (light gray dot). The horizontal dashed lines roughly indicate the 95% 
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for 0.2 < z < 0.45 and 0.45 < z < 0.7 from Noes ke et alj <2007J). The green 
dotted curve indicates the "quenching mass limit" as presented in Bu ndv etal] 
1 2007). (A color version of this figure is available in the online journal.) 



TABLE 2 

Detection Rates for Na i D Outflows Selected by Lfr 



Criterion 




^outflow 


Detection Rate 


Within MIPS coverage 


169 


24 


0.14 ±0.03 


log(L IR /L )> 11 


21 


8 


0.38 ±0.11 


log(L IR /L )< 11 


64 


5 


0.08 ±0.03 


No 24 fi detection 


84 


11 


0.13 ±0.04 



NOTE. — The conventional cut for LIRGs is log(LrR/Z,Q) > 
11. Only the objects with high-S/N Nal D velocities are included 
for the calcula tion of detection rates; i.e., the objects in the low- 
S/N sample (§ |2~4~fl are treated as if they are not detected. The 
uncertainties are estimated from binomial statistics. 

Some insights into these effects are gained from Fig. [6] Ap- 
parently, the high-S/N sample (§ 12. 4. It does not come close 
to the completeness limit in stellar mass of the parent AEGIS 
survey, except at the lowest redshifts. The lower sampling 
rate of high-mass galaxies (e.g., log (M*/M Q ) > 11) at low 
redshift may be a result of the smaller survey volume; that 
is, lu minous early-ty pe galaxies tend to be highly clustered 
(e.g. JCoil et al.f 2008). and a smaller number of overdense re- 
gions with bright ellipticals fall in the field of view toward 
lower redshift. The effect of the changing physical aperture 
size, from 2 kpc to 9.5 kpc corresponding to the 1" slit, on 
the detectabil ity of outflow is difficult to assess. The work by 
Martin ( 2006) on the spatially resolved Na I D outflows in lo- 
cal ULIRGs indicates that the extended blueshifted interstel- 
lar absorption is found over the scale of > 15 kpc, in which 
case a significant column should remain available within the 
regions covered by the slit over the redshift range covered in 
this study. 

3. HOST GALAXIES OF OUTFLOWS 

3.1. Trends with Star Formation 

Since the recent advance in the detailed knowledge of 
starburst-driven galactic winds have been gained through 
the studies of local LI R Gs (e.g.. IHeckman et alj 120001 : 
Rroke et alJl2002L 1200535 |Martin 200§, the NaiD veloc- 
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Fig. 7. — NaiD velocity as a fun ction of infrared luminosity for the 
objects in the high-S/N sample (j) |2.4.U . The figure includes only the objects 
that are covered by the Spitzer/MlPS observations and are detected above the 
flux limit (524 ~ Jy) to yield the measurements of Lir. The symbols are 
as in Fig. |2] (A color version of this figure is available in the online journal.) 

ity measurements as a function of infrared luminosities Lir 
naturally provide a starting point for comparison. Since for 
star-forming galaxies the infrared is dominated by the ther- 
mal dust emission of reprocessed starlight from hot, young 
massive stars, a tight correlation e xists between L ir and 
SFR in dusty star-forming galaxies (Kennicutt 1998). For 
a subset of our sample, the far-infrared photo metry from 
5p;tee r/Multiband Imaging Photometer (MIPS; Rie ke et alj 
|2004|) i s used to derive the tot al infrared lum inosity Lir, fol- 
lowing iLe Floch etal] (120051) and using the IChary & Elbazl 
(2001) spectral energy distribution (SED) templates. 

In Fig. [7] a clear tendency is observed for high Lir objects 
to host outflows. Indeed, a majority of outflows are found 
in LIRGs (Lir > 10 11 L©) in the figure. Table|2]presents the 
outflow detection rates for the subsampling based on Lir. The 
outflow detection rate of38%±ll% for LIRGs is similar to 
those reported by low -z surveys of infrare d-selected galax- 
ie s; e.g., 42% ± 8% inlRupke et alj d2005bl) and 32% ± 12% 
in IHeckman et al . (2000). Hence we confirm the results re- 
ported by others that the detection rate of outflow in infrared- 
selected galaxies correlates well with their infrared luminos- 
ity. We note, however, that the exact values for the de- 
tection rates, especially among the galaxies with lower Lir, 
may not be robust against selection effects and incomplete- 
ness, since our primary selection cut is made by the strength 
of continuum around Nal D (§ 12.31 ). There are no ULIRGs 
[log(LiR/L Q ) > 12] in our sample, which is reasonable given 
the expected number of ULIRGs in the survey volume is very 
small (< 10), assuming r edshift-dependent l uminosity func- 
tions of infrared sources dLe Floc'h et alj|2005l) and the com- 
pleteness of the DEEP2 survey. 

In Fig. [8] the Na I D velocity as a function of the total SFR 
(i.e., the sum of SFR s derived from infr ared SED and optical 
emission lines; see Noes ke et al.l ((2007a) for details) is shown. 
The similar dependence of v(Na I D) on Lir and the total SFR 
is expected due to the tight correlation between Lir and the 
total SFR, which is often dominated by the infrared contri- 
bution. Again, a majority of outflows are seen in galaxies 
with SFR > 20 M Q yr" 1 , which correspond to the amount of 
SFR expected in LIRGs. As a caveat, it should be mentioned 
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FIG. 9. — Specific SFR as a function o f ste llar mass for blue-cloud (top) 
and red-sequence (bottom) galaxies; see j) l3.2l for the color cut scheme. The 
symbols are as in Fig.|I] The fi gure is comparable to Fig. 1 of Noe ske et alj 
1 2007b), except that red AGN/LINER candidates are not removed. An SFR 
estimate from optical line emission suffers significantly from AGN/LINER 
contamination I Yan et al. 2006; Weiner et al. 2007); therefore, their specific 
SFRs are likely overestimates here, especially for the red-sequence galaxies. 
(A color version of this figure is available in the online journal. ) 

that some of the infrared luminosity might be from (obscured) 
infrared-bright AGN, and not from star formation. We lack 
the diagnostics to distinguish star formation from AGN, but 
the simple relation between Lir and SFR may not hold for the 
cases in which AGN contributions to the infrared luminosities 
are high. Local studies of LIRGs suggest that strong AGN 
contamination may be s mall, though low- luminosity AGN 
could be common (e.g.. lVeilleux et al.lll995l) . 

Since t he SFR scales almost p roportionally with galaxy 
mass (e.g., Brinchma nn et al.l2004l) . its inference with respect 
to the specific amount of star-forming activity m ay be mis- 
leadin g. The birthrate parameter, b = SFR/ (SFR) dKennicuttl 
fl998h . is a better indicator which takes into account the star 
formation history yet is difficult to estimate due to its depen- 



dence on the timescale over which the galaxy has been form- 
ing stars as well as on the fraction of gas recycled in the past. 
It is still desirable to remove the first-order effect of galaxy 
mass, since our sample is drawn over the redshift interval 
where downsizing affects star formation particularly strongly. 
Specific SFR, defined as SFR/M* (i.e., SFR per unit stellar 
mass) is easier to compute and a reasonable proxy for birthrate 
parameter. Fig.|9]shows how the specific SFRs are distributed 
as a function of galaxy stellar mass M», obtained from fitting 
SEDs to optical and near-infrared photometry (Bun dvet al.l 
2006). For the red-sequence galaxies (Fig.QtjJl, the specific 
SFRs may be overest i mated due to AGN /LINER contamina- 
tion (lYan etal.ir2006t IWeiner et aI1l2007l) . Focusing on the 
blue-cloud galaxies, an interesting feature is that, at a fixed 
mass, a majority of the objects with outflows are seen at the 
upper envelope of the distribution of specific SFRs. The ana l- 
ysis of the parent AEGIS sample by Noes keet al.l (l2007allbh 
indicates that, in the redshift interval 0.11 < z < 0.55 from 
which our sample is drawn, M» ~ 10 11 M© is roughly where 
star formation is se en to be quenched in the most massive 
galaxies. In the full lNoeske et al.l (l2007bl) sample, the distri- 
bution of star-forming galaxies is fairly tight; the loci in Fig. [9] 
suggest that host galaxies of outflows are undergoing an en- 
hanced episode of star formation, relative to the star-forming 
galaxies without outflows at the same epoch. 

It is emphasized that star formation may not be the only 
source of optical/infrared emission in the presence of an AGN. 
We see that quite a few red-sequence galaxies with inflows 
show up in the figures involving SFRs presented in this sec- 
tion. As discussed in the next section, these objects show 
little evidence of ongoing star formation in UV diagnos- 
tics, and most are either quiescent or show Seyfert/LINER- 
like excitation in their [Nil] /Ha emission line ratios (§ 14. 31 . 
Hence inflows in general are more likely to be associated with 
AGN/LINERs than star formation. 

3.2. Color Magnitude Diagram 

Although optical photometry only tells an incomplete his- 
tory, dividing galaxies into blue and red populations by op- 
tical colors still is very useful for capturing the essence of 
galaxy evolution with a bird's-eye view, since the bimodal- 
ity of galaxy color distribution is among the most prominent 
features that persist over wide ranges of parameters, such 
as galaxy mass, luminosity, environment, and redshift (e.g. , 
Baidr v~etalll2004l iBalogh et all 120041: IWillmeret all l2006t 
Coop eret al.1 120071) . Over the past few decades, numerous 
studies of local galaxies have shown that the "red sequence" 
is generally populated by red, dead, passively evolving galax- 
ies with old stellar populations, while the "blue cloud" is 
populated by blue, actively star-forming galaxies. These two 
densely populated regions are divided by a sparsel y populated 
area, sometimes called the "green valley"; see iFaber et all 
d2007l) and references therein for a comprehensive account. 

Fig.lTOlpresents the rest-frame (U-B,Mb) color-magnitude 
diagram, in which the host galaxies of outflows are seen in 
relation to the red sequence and the blue cloud. The cata- 
log of rest-frame pho tometry was constructed as described in 
Willm er" ;t al.l ((2006). 1 1 Given that the low-z red sequence 
galaxies are characterized as inactive and quiescent, that a 
majority of outflows are found in red galaxies comes as a sur- 
prise; among 32 outflows, 21 of them are found on the red 



11 The optical photometry are on the Vega system. See Will mer et al] 
1 2006 ) for the conversion procedure between Vega and AB magnitudes. 
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FIG. 10. — Rest-frame (U — B)-Mb color-magnitude diagram. The pho- 
tometry, on the Vega system, is correc ted for Galactic ext inction but not for 
the internal extinction of galaxies; see Willmer et al. 1 2006) for detail. The 
symbols are as in Fig. [2] The dashed line indicates the U -B color cut used 
to divide red-sequence and blue-cloud galaxies; see Faberetal. 12007). [See 
the electronic edition of the paper for a color version of this figure.] 



sequence. Non-star-forming, quiescent galaxies are red in the 
rest-frame optical due to their intrinsic SEDs; the rest-frame 
U -B in particular brackets the 4000 A feature, which is sen- 
sitive to the "break" caused by the Balmer limit, as well as the 
high metal opacity in the cool stellar atmosphere of old stars. 
Nonetheless, high dust reddening can push some star-forming 
galaxi es into the red sequence, especially at higher redshift 
(e.g., IWeiner et alj l200l iBell et al.1 120051) . The rest-frame 
photometry in Fig. [I0]are not corrected for the internal extinc- 
tion of galaxies, so the "contamination" from dust-reddened 
star-forming galaxies can happen though expected to be small 
at z < 0.5. We also reiterate that, due to the fixed contin- 
uum S/N cut, our sample is naturally biased for red-sequence 
galaxies, which tend to have high surface brightness (§ !2.31 l. 
This selection bias due to the continuum strength can be in- 
ferred from the difference in the outflow detection rates for red 
and blue subpopulations. Among the full sample (/V = 431), 
the outflow detection rates are 1 1 /205 = 0.05 ± 0.02 and 
21/226 = 0.09 + 0.02 for blue and red objects, while they 
are 1 1 /5 1 = 0.22 ± 0.06 and 21 /154 = 0. 14 ± 0.03 among the 
high-S/N sample (N = 205). 

Fig. QT] shows the color-magnitude diagram in which the 
objects are denoted by their infrared luminosity and offers 
evidence that star-forming galaxies do reside in the red se- 
quence, though relatively small in fraction. The sample is 
divided at z = 0.35 into low- and high-redshift subsamples, 
which makes clear that the infrared-luminous galaxy fraction 
is much greater at higher redshift, consistent with lookback- 
time studies of star-forming, infrared-luminous galaxies (e.; 
IBell et al.ll200l iLe Floc'h ei^ 12001 IKfoeske et al.1 12007a 
Outflows are mostly observed in z > 0.35 objects. A com- 
parison of the top and bottom panels shows quite strikingly 
that the z > 0.35 red sequence has a substantial number of 
infrared-luminous galaxies, some of which are LIRGs, which 
is absent in the z < 0.35 red sequence. The red-sequence out- 
flows, however, are also found in the objects without signifi- 
cant 24 fim flux. 

While the narrow spectral baseline of DEEP2 makes the 
Balmer decrement unavailable, the UV spectral slopes (3 from 
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FIG. 11. — Rest-frame (U-B)-Mb color-magnitude diagram for objects at 
z > 0.35 (top) and z < 0.35 (bottom). The gray dashed lines divide red and 
blue galaxies (Fig.[l0j. The symbol indicates the infrared luminosity Lir: 
below flux limit (dot), log(Li R /L Q ) < 10.5 (cross), 10.5 < log(L IR /L Q ) < 
11 (open circle), and log(LiR/LQ) > 11 (open square). The color scheme 
for the symbols is similar to Fig.[2]and indicates Nal D kinematics: outflow 
(blue), inflow (red), systemic (gray), and low S/N (light gray). There are 
roughly equal number of objects above and below the redshift cut at z = 0.35. 
[See the electronic edition of the paper for a color version of this figure.] 
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FIG. 12. — Rest-frame U—B color as a function of the UV spectral slope /3 
for a subset of objects with GALEX photometry. The (i follows the definition 
given by Seibert et al. (2005), f\ oc \P; star-forming galaxies have smaller 
values of ft due to their rising continua toward shorter wavelengths. The 
symbols are as in Fig. [2] A majority of blue-sequence galaxies are classified 
as low-S/N; they are bright in UV but fainter in the visible, making it difficult 
to obtain their high-S/N spectra. (A color version of this figure is available in 
the online journal.) 

the Galaxy Evolution Explorer (GALEX) photometry are 
measured for a subset of our sample. 12 The UV spectral slope 
is generally a good measure of dust reddening, since it is a 
direct measure of continuum slope in ultraviolet, where dust 
optical depth is particularly high. Fig. [121 shows that some 
red galaxies with outflows have j3 < 0, comparable to typ- 
ical blue star-forming galaxies. None of these with 24 /im 
detection, however, are LIRGs. Another look into their iden- 

12 The images were processed using ver. 4.1 of GALEX pipeline; see 
http://www.galex.caltech.edu 
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Fig. 13. — Top: rest-frame U —B color as a function of rest-frame 
NUV— Rab color. The horizontal dashed line al U — B = 0.15 roughly di- 
vides red-sequence and blue-cloud galaxies (Fig. 1101 . The vertical dashed 
line at NUV-Rab =5.4 follows the delineat ion of star-f orming and non-star- 
forming galaxies employed by Schawinski et al. ( 2006), who used NUV — r; 
r is roughly comparable to Rab ■ The symbols are as in Fig.|2] Bottom: H/3 
spectral line index as a function of rest-frame NUV — Rab color for a sub- 
set of objects with the H/3 coverage. The strength of H/3 absorption [i.e., 
Wb(H/3) > A] is sensiti ve t o a starburst occurred within ~ 10 7 Myr-2 Gyr 
ago (see the inset of Fig. [141 . A negative Wo(H/3) indicates the presence of 
emission at H/3, i.e., Wq is an estimate of the sum of emission and absorption 
equivalent widths (§ 12.2. 1> . (A color version of this figure is available in the 
online journal.) 

tity is given by the top panel of Fig. Qj] where the distribu- 
tion of galaxies in their rest-frame U—B color is compared 
to that of the rest-frame NUV-/?ab color. The UV-optical 
color is much more sensitive to a small fraction of young 
stellar populations and in turn the current as well as recent 
star-forming activity than optical colors alone. The particular 
NUV-^ab cut delineating star form ation from no star forma- 
tion is from lSchawinski et al. ( 2006), who used NUV-r col- 
ors; SDSS r is close to CFHT Rab- In this empirical color cut, 
early-type galaxies with the strongest UV upturn observed lo- 
cally should not contaminate the star-forming classification; 
this makes inevitable that some objects in the "no star for- 
mat ion" region may actually be star forming. Based on the 
cut. lSchawinski et all found that ~ 30% of visually classified 
z < 0.1 early-type galaxies brighter than M r = -21.5 showed 
signs of recent star formation. Strikingly, most optically red 
galaxies with outflows have NUV-/?ab < 5.4, suggestive of 
recent star formation. 

These results give rise to an interpretation that the ar- 
rival on the red sequence of the galaxies with outflows hap- 
pened only recently; i.e., red sequence outflows are found 
predominantly in post-star-forming galaxies with detectable 
residual star formation or dusty star-forming galaxies with 
high infrared emission. The redness in their visible colors 
might also arise partly from the presence of dust; the strong 
correlation between reddening E(B-V) and the equivalent 
width of low-ionization absorption lines are often reported 
fe-g JArmus et al.lll989b IVeiUeux et allll995t iHeckman et all 
2000). Poststarburst galaxies, identified spectroscopically 
by their strong B aim er absorption, are generally found to 
be dusty as well (e.g., iPogg ianti & Wu 2000; Pog gianti et al.1 
l200lHBalogh et alJfeoolUSato & Martini 120061) . The appear- 
ance of the post-star-forming phase also fits well with the 
scenario that it is observed when star-forming galaxies in the 
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FIG. 14. — Nal D velocity as a function of H/3 spectral line index. The 
symbols are as in Fig. [2] Since the H/3 index measures the sum of absorp- 
tion and emission lines, the H/3 index becomes a lower limit to the strength 
of H/3 absorption equivalent width, if emission filling is significant. In turn, 
Wb(H/3) < A indicat es tha t the emission line flux is greater than that of 
absorption line. See § 12.2. 1 1 for detail on the definition of the spectral line 
index. Inset: the stellar spectral line index at Ha (red dashed line) and H/3 
(blue solid line) as a function of stellar age s. The indices are measured from 
a single stellar population model from Delgado et al. (2005). A stellar pop- 
ulation with the age > 10 Gyr has W (H/3) fa 2 A. An H/3 index > 3 A is 
expected for a stellar population during a poststarburst phase from ^10 Myr 
up to ~ 2 Gyr. (A color version of this figure is available in the online jour- 
nal.) 

blue cloud make the transition to the red-sequence after some 
mechanism triggered an enhanced episode of star formation, 
which then gets shut off. The NUV-^ab color alone, how- 
ever, does not rule out low-level star formation in these galax- 
ies, since whether a galaxy goes through poststarburst de- 
pends on the fraction of mass that has formed in the most 
recent star formation event as well as its timescale. It is also 
worth noting the clear separation between outflows and in- 
flows in NUV-^ab; most inflows are observed in red galaxies 
in which little or no star formation is detected. 

3.3. Evidence for Poststarburst 

The Balmer absorption lines are sensitive to the age of 
the underlying stellar population and become prominent in 
the spectrum in which A to early-F stars, living up to ~ 
1.5 Gyr, contribute significantly. This feature has been ex- 
ploited to find the signature of poststarbursts in distant galax- 
ies, often using HS and/or H7 absor ption due to less con- 
tamination from emis s ion filling (e.g. , [Dressier & Gun nl 19831 : 
IZabludoff et aH[l996l iDressler et alJl2(M Ig^ 120071) . The 
limited spectral baseline makes only low-lying Balmer lines 
available for our sample. Fig.[T4lshows the Nal D velocity as 
a function of the H/3 spectral line index Wo(H/3). A spectral 
line index is an estimate of the sum of emission and absorp- 
tion equivalent widths (§ 12.2.11 ). In the figure, outflows are 
seen in the following two classes of objects. One class is those 
with Wb(H/3) < A, meaning that H/3 is seen in emission, so 
they are the outflows seen in star-forming galaxies in the blue 
cloud (Fig. [Toll. Another is those with Wq(U/3) > 3 A, slightly 
offset from a concentration of objects with systemic Na I D 
velocity around Wo(H/3) «2 A. To put this in perspective, 
the inset of Fig.[T4lshows the evolution of stellar H/3 absorp- 
tion index over a range of stellar ages. A stellar population 
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with the age > 10 Gyr has its W»(H/3) asympototing to m 2 A. 
The Hf3 absorption would be seen at Wb(H/3) > 3 A during 
the poststarburst phase from 10 Myr to ~ 2 Gyr. In prac- 
tice, a measured H[3 index from a galaxy spectrum includes 
some nebular emission line flux originating from the H II re- 
gions surrounding hot, young stars in the presence of resid- 
ual star formation, pushing a Balmer absorption index to a 
smaller, less positive value. Therefore, Wb(H/3) > 3 A objects 
are likely to be in the poststarburst phase. A large number of 
systemic Na I D velocity objects around Wq(H(3) « 2 A is also 
consistent with old, quiescent galaxies not hosting outflows. 

The bottom panel of Fig. [l3]shows the relation between the 
NUV-/?ab color and the H/3 spectral line index for the sub- 
set of sample with both measurements. Although the small 
statistics make interpretation difficult, we can see that the 
red-sequence objects with outflows appear to have a slightly 
higher Wb(H/3) in general compared to those without. A few 
blue-sequence objects having NUV-/?ab < 5.4, both with 
and without outflows, also move to the Wo(H/3) > A regime. 
If these objects host residual star formation, inferred from 
their blue NUV-/?ab color, their emission fluxes may fill H/3 
absorption, making the observed Wo(H/3) w 2-3 A only lower 
limits to their absorption equivalent strength. This would 
make the case stronger for the poststarburst identity. 

The combination of selection cuts from our sample yields 
a small number of fully overlapping objects across a variety 
of measurements and makes it difficult to reach generalizing 
conclusions with statistical rigor. Nevertheless we do find 
overwhelming consistency in the evidence that our LIRG-type 
outflows are mostly seen either in starburst or poststarburst 
objects. This naturally fits into the currently favored scenario 
of galaxy evolution between blue-cloud and red-sequence 
galaxies, in which some mechanism (e.g., merger) triggers 
a starburst in a blue galaxy which then gets "quenched" by 
some feedback mechanism, such as by an AGN or super- 
novae, by the time the galaxy joins the red sequence. The 
outflows may be the result of such feedback process observed 
in "transition" objects. 

However, we must also note that stronger Balmer absorp- 
tion only indicates that there was certainly a detectable en- 
hancement, relative to the present, of star formation in the 
recent past; whether or not an individual case makes the cri- 
teria for a conventional starbusrt is admittedly unclear, espe- 
cially given our crude diagnostics. There are plausible ways 
for galaxies to quench star formation without going through 
a starburst phase, and such mechanisms could generate out- 
flows and detectable enhancement of Balmer absorption, if 
quenching occurs quickly. This possibility can be explored 
further only through better diagnostics. 13 

3.4. Host Morphology 

Fig- US shows a subset of galaxies with NalD measure- 
ments for which we have the quan titative morpholo gy from 
the HSTI 'ACS imaging analysis by lLotz et aT] ( 120081) . A set 
of new nonparametric morphology measures, Gini coefficient 

13 With our sample, no robust way exists for directly connecting poststar- 
burst to the kind of star formation detected in NUV. The advantage of NUV 
diagnostics here is its sensitivity to low-level star formation that normally 
eludes detection in the visible, due to much larger light contribution from old 
stars to the visible region of galaxy spectrum, and may not tell us much about 
the timescale over which such low-level star formation has been in existence. 
Therefore the star formation seen in NUV could be "residual" star formation 
from a recent starburst/star-forming event or it could just be a small amount 
of continuous star formation. 
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FIG. 15. — Quantitative morphology measure (Gini coefficient G and 
the second-order moment of the brightest 2 0% of the tot al galaxy flux M20) 
from the analysis of HSTI ACS images by Lotz et al. 1 200S§). The morpholog- 
ical cu ts (dashed gray lines) to the "traditional" classes are from lLotz et all 
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ple. The symbols are as in Fig. [2] (A color version of this figure is available 
in the online journal.) 
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circle), late-type (Sb-Ir; square), and merger-candidates (star). The filled 
symbols indicate host galaxies of outflows. The symbol colors are as in 
Fig. [2] Only the objects from the high-S/N Nal D velocity sample (§ 12.4.11 
with NUV measurements are plotted. (A color version of this figure is avail- 
able in the online journal. ) 

G which measures the distribution of flux among pixels over 
a galaxy and the second-order moment of the brightest 20% 
of the total galaxy flux M20, is calibr ated visually to th e tradi- 
tional morphological classifications dLotz et a l. 2004): early- 
type (E/SO/Sa), late-type (Sb-Ir), and merger candidates. Our 
survey is more sensitive to luminous, high surface brightness 
galaxies (§ 12.31 ). and Fig.[l5lshows that the sampling is biased 
against objects w ith late-type (i.e., lower surface brightness) 
morphology; see lLotz et ail (120081) for the analysis of com- 
plete samples drawn from the parent AEGIS survey, which 
shows that the many objects with low surface brightness, late- 
type morphology, which would appear toward the bottom-left 
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corner, did not make our selection cut. 

The transition between a LIR G to a ULIRG is physically 
plausible in the merger sequence (Sanders 2004), so the high 
outflow detection rate (3 /6 = 0.5 ± 0.2) in merger candidates 
is not surprisin g. Outflows are detected in almost all lo- 
cal ULIRGs dMartin| HH Hop EuBjEeLaD |2005b). Al- 
though all three merger candidates with outflows are de- 
tected at 24 /im, only one is a LIRG. Since the dynamical 
timescale of merger, ~ 1 Gyr, exceeds that of gas consump- 
tion, ^100 Myr, it is likely that we miss the (U)LIRG phase 
due to its short duty cycle. It should be noted, however, the 
nature of the association between mergers/interactions and 
(U)LIRGs is not as well established for objects at z ~ 0.5. 

A majority (9/14 = 0.64 ±0.13) of outflows are seen in 
the objects with early-type morphology (E/SO/Sa). An in- 
sight on the order of morphological transformation may be 
given by Fig. [16] which plots the distribution of morphology 
in terms of the rest-frame ({/-B,NUV-/?ab) color-color di- 
agram. Consistent with the common knowledge, most late- 
type galaxies are in the blue cloud, while most early-type 
galaxies are in the red sequence. The distribution of outflow 
hosts roughly follows that of the parent sample. Although the 
small sample size does not allow a statistically robust con- 
clusion, outflows appear to be seen more preferentially where 
merger candidates are also seen. Previously we find evidence 
that some outflows are seen in poststarburst objects (§ !3.31 l. 
The early-type morphology, coupled with the signs of inter- 
action in a few of the red-sequence outflows, is consistent 
with the existing studies of poststarburst galaxies as transi- 
tion objects in the morphology sequence, which conclude that 
they are spheroidal, often sh owing signs of interaction (e.g., 
lYang et al.ll2004t lGotoll2005h . 

We find only two outflows with late-type morphology, at 
Z > 0.35. Fig. [T5l makes clear that our selection is highly bi- 
ased toward high surface brightness galaxies; we have missed 
a large population of low surface brightness star-forming 
galaxies that are dim in the rest-frame visible continuum, 
with which our primary selection cut was made. These out- 
flows with late-type morphology are detected significantly in 
the infrared, having log(LiR/L Q ) w 11.1 and 10.4. Unlike 
ULIRGs, the morphology of LIRGs appears to be a mixed 
bag. At z ~ 0.7, however , more than half of their LIRGs 
have disk morphology (e.g.. lBelfet al. 2005; Mel bourne et al.1 
2005), suggesting that the abundant LIRG population simply 
reflects an elevated level of star formation in normal galaxies 
at an earlier lookback time. All our objects are at lower red- 
shifts, and it is unclear if we are observing disky LIRGs of the 
type seen at higher redshifts. 

4. DISCUSSION 

We have seen that the detection rate of LIRG-like outflows 
is a st rong function of infrared luminosity or star formation 
rate (§ 13. 11 1. The distribution of specific star formation rates 
at a fixed stellar mass (Pig. |9]» might imply that the frequency 
of outflows may be higher for the objects in which ongoing 
star-forming activity is enhanced relative to the past average 
(i.e., high birthrate parameter). The distribution in the opti- 
cal color magnitude diagram shows that outflows are abun- 
dant both in the blue cloud and the red sequence (§ I3.21 >. The 
NUV photometry reveals that the red-sequence outflows are 
hosted by the objects that have gone through recent star for- 
mation (Fig. [T3Tl. There is evidence that some of the red- 
sequence outflows are in dusty (Fig.[T2l. star-forming galax- 
ies (Fig. |TT}. Yet some may be poststarburst galaxies (§ !3.31 l 



just arriving the red sequence, which is consistent with their 
early-type morphology (§ 13. 41 . We now discuss a few fine 
points that need attention, and describe how our finding of 
Na I D outflows might fit into the current understanding of 
galaxy formation and evolution. 

4.1. Redshift Evolution of Outflows 

Fig.|6]shows the detection of outflows as a function of red- 
shift in our sample. Over the redshift range, star-forming 
properties of galaxies are expected to change and so do the 
detection rates of outflows. In the figure, no M* rj 10 n M Q 
galaxies at z ~ 0.2 show signs of an outflow, yet a signifi- 
cant fraction of galaxies with the similar mass at z ~ 0.5 do. 
This is consistent with the known trend of the star formation 
history of galaxies an d the increasing trend of ou tflow detec- 
tion rate with SFR (§EB. iNoeske et all (l2007al) shows that 
star-forming galaxies with M* ss 1O U M0 are rare at z < 0.45 
but are abundant and most are LIRGs at z > 0.45. Fig. [6] 
shows that the outflows detected in our sample are typically 
seen in high-M* galaxies at higher z with an SFR expected 
in a LIRG (> 20 M* yr" 1 ). Due to low S/N, however, we 
cannot say whether or not lower-M, galaxies at z ~ 0.5 host 
outflows. Nevertheless, at the high M* end, we clearly ob- 
serve that the "downsizing" effect extends to the detection 
rate of LIRG-like outflows, i.e., the typical mass of galax- 
ies that host outflows move to a lower mass toward lower 
redshift, much like the trend seen in cosmic star formation 
history of galaxies (e . g..lMadau et alJll996tlLillv et al. ll 19961; 
Le Floc'h et alJl2005t iFaber et alJl2007t INoeske et al] l2007at 
Cooper et al . 2008) . iBundv et a l. (2006) qunatify the evolving 
trend in terms of the "quenching mass limit," which is shown 
in Fig. [6] for a crude comparison. An important connection 
is likely to exist between quenching events and gaseous feed- 
back, so the redshift evolution of outflows needs to be ex- 
plored in view of host galaxy properties. 

4.2. Merger-Triggered Activities and Outflows 

Merger-driven galaxy evolution has been a prominent 
paradigm over the past decades, and the current success of 
the ACDM theory in describing the hierarchical structure for- 
mation certainly points to its significant roles in galaxy-scale 
phenomena. The relevance of mergers in shaping the cos- 
mic star formation history of galaxies, however, is under in- 
tense scrutiny. The difficulty arises from the fact that di- 
rect identification of mergers and quantifying their frequency 
from the observation of high-redshift galaxies remains very 
challenging, due to a number of factors including surface 
brightness dimming and shifting passbands. Galaxy mergers 
of star-forming disk galaxies have attracted much attention, 
because they are a mechanism widely known to cause en- 
hanced star formation, after which merger remnants dynam- 
ically relax into spheroids (e.g., Mihos & Hernquist 1996; 
ICoxetai1 l2006). Along with the expected increase of the 
merger rat e with the lookback time from numerical simula- 
tions (e.g., Gottltiber et al. 2001), mergers have been invoked 
as a mechanism responsible for the declining tren d of the co- 
moving star formation d ensity since z ~ 2 (e.g., iLillv et aTl 
Il996t iMadau et afl fT996). although exactly how much merg- 
ers contribute to the tr end is still debated (e.g., Brid ge et all 
l2007llLotzet"aT]l2008h . Galaxy merger is appealing also for 
it affects star-forming as well as morphological properties of 
galaxies in ways that naturally explain the observed transi- 
tion of young (i.e., blue and disky) into old (i.e., red and 
spheroidal) galaxy populations. In principle, a variety of such 
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transient phenomena as starburst and quasars can also be in- 
tegrated int o a coherent picture of merger-driven galaxy e vo- 
lution (e.g. JSanders & Mirabellll996tlHopkins et al.ll20 06). 

The correlation of the outflow detection rate with the de- 
gree of elevation in star-forming activity certainly indicates 
that starbursts are an important component of galactic super- 
wind phenomenon. However, outflow is exciting not because 
it is associated with such a spectacular starburst event, but 
because it may carry away from the host galaxy the bulk of 
"fuel" for further star formation, which may contribute to the 
"quenching" of star formation. In fact, ULIRGs (i.e., gas-rich 
mergers) show dynamical ev idence for spheroids in forma- 
tion (e.g., iGenzel et~aT1l2001l) . and the star formation history 
of post starburst galaxy is consistent with the (LHLIRG origin 
(e,L iPoggianti & Wull2000t iBekki et alJl200U iKavirai et alj 
120071) . Our detection of NalD outflows in poststarburst 
galaxies strongly suggests that feedback mechanism affects 
the kinematics of cool interstellar gas well after the most in- 
tense phase of star formation. The fact that quite a few red- 
sequence galaxies host winds may not be surprising yet still 
a striking result. A vast majority of absorption-line studies 
of outflows in literature has been on vi gorously star-formin 
systems. Using a plot similar to Fig. fTT lRupke et all (12005 
showed that most outflows in infrared-selected galaxies were 
detected in the loci of low Mg I b and high Na I D absorption 
indices, i.e., young galaxies with high interstellar Nal D col- 
umn (their Fig. 8). In contrast, the red-sequence galaxies with 
outflows in our sample can have a high Mg I b absorption in- 
dex, reflecting the presence of intermediate to old stellar pop- 
ulation, and are not clearly distinct from other red-sequence 
galaxies in an optical color magnitude diagram (Fig.lTOt. The 
apparent connection between outflows and high M* poststar- 
burst galaxies provides circumstantial evidence that outflows 
of cool gas contributes to or is a consequence of more effec- 
tive quenching of star formation. 

How the progenitors of spheroidal galaxies in the local 
universe evolve into their current state remains a topic un- 
der vigorous investigation. While their stellar contents sug- 
gest that massive spheroids have been passively evolving 
and that their stellar mass changes little sin ce z ~ 1 (e.g., 
iBrinchmann & Ellis! 120001 iBundy et al.ll2005l) . the evolution 
of the luminosity function suggests that the number den- 
sity of luminous red galaxies has increased by a factor of 
~ 2 oyer the same peri od (e.g., iBell et al.1l2004t iBrown et"ai1 
l2007l;lFaber et al.ll2007l) . Over the similar redshift range, star- 
forming galaxies have a relatively smal l spread in SFRs at a 
fixed mass in the blue sequence (e.g.. iNoeske et al] l2007a). 
leading to a paucity of objects presumably in transition be- 
tween the blue and red sequences. Furtherm ore, less massive 
spheroids haye you nger stellar contents (e.g jTreu et al.l2005t 
iKaviraj eTa l. 2008) and the characteristic galaxy mass above 
which the star formation in gala xies quenched evo lves over 
redshift in a downsizing fashion (Bundv et al. 2006), indicat- 
ing that catastrophic transition events occur at a progressively 
lower mass toward lower redshift. The exact rate of transi- 
tion is very difficult to estimate, yet indirect argum ents favor 
rare and/or fast mechanism(s) (e.g.. lBlantonll2006l) . Gas-rich 
mergers perhap s contribute to som e but not all of these blue- 
red transitions (Bundv et al. 2007). The qualitatively similar 
downsizing trend in star-forming galaxies and outflow hosts 
over z < 0.5 (§ 14. j} , however, may imply that the mechanism 
responsible for downsizing of star formation may also accom- 
pany outflows. A rigorous conclusion must await the quanti- 
tative analysis of the sample which suffers less from small 



number statistics, selection effects, and incompleteness. 

We must also note that it is not clear that merger per se 
is a necessary precursor for outflows. The direct morpho- 
logical evidence for interaction is not very strong in our out- 
flows (a majority of outflows are in early-type galaxies; § I3.41 >. 
A circumstantial evidence is provided only through an indi- 
rect argument that a spheroidal formation follows a merger- 
induced starburst with a poststarburst signature. Since the 
timescale for mergers at high redshift to remain identifiable 
is shorter than a poststarburst phase, the scarcity of direct evi- 
dence may not immediately discount the importance of merg- 
ers. However, it is possible that low-level star formation in 
early-type galaxies could drive outflows, perhaps via the ac- 
cretion of gas-rich satellites or minor mergers. Early-type 
galaxies also have reservoirs of hot gas, which could provide 
fuel for star formation via condensation. It has been suggested 
that the early-type galaxy population itself shows bimodality 
in their UV-visible color distribution, reflecting the richness in 
their star formation histories; low-level star forma tion appears 
to con tinue in less-massive early-type galaxies (IKaviraj et al.1 
2008). The origin of outflows in the red galaxies may not be 
as simple as a quenching event followed by passive evolution. 

Furthermore, the mounting evidence now shows that a ma- 
jority of z ~ 1 LIRGs are normal disk galaxies whos e elevated 
star formation is not caused by interaction (e.g., IBell et all 
l2005t lLotz et aO2008l) . The gradual decline of star-forming 
efficiency in their disks may be responsible for much of the 
global trend seen in comoving SFR density. From our study 
alone, whether or not these disky LIRGs at high-redshift host 
outflows is not clear; the nature of outflows in these objects 
may be different from merger-induced ones. 

It would be interesting to see how the presence of out- 
flows at z < 1 contribute to the fate of z ~ 1 LIRGs down 
to z ~ 0, which may either stay but fade gradually within the 
blue sequence or go through rapid quenching of star formation 
to migrate to the red sequence. Given the strong evolution 
of galaxy properties in general (§ 14. Il l, our knowledge from 
the local study of LIRG-like outflows might not be relevant 
for high-z, disky LIRGs. On the other hand, if the outflow 
strength correlates with some parameter such as the presence 
of nuclear activity or their degree of interaction in relation to 
their morph ology among h igh-z LIRGs [as observed in local 
ULIRGs by Martin (2005)], deprivation of (cool) gas via su- 
perwind may be important in transforming star-forming disks 
into quiescent spheroids in catastrophic events at those red- 
shifts. The detection rate of outflows in view of host galaxy 
morphology at z < 1 may provide some insight on the physical 
mechanism that maintains the bimodality in galaxy popula- 
tion since z ~ 1. A comparison of mass outflow rates in z ~ 1 
LIRGs with different morphology would also make an inter- 
esting exercise to put some constraint on the role of mass load- 
ing in the evolutionary histories of morphology and star for- 
mation. Fortunately, several useful UV resonance lines at dif- 
ferent ioni zation states shift in to visible window for z ~ 1 ob- 
jects (e.g., Weiner et al. 2009), so future surveys are in a bet- 
ter position to con strain mass loading from these lines (e.g., 
iMurray et al.ll2007l) . 

4.3. Star Formation Versus AGN 

It is usually assumed that superwinds are driven by ther- 
malized energy output from supernovae. In principle, AGNs 
offer much larger repository of energy than supernovae, and 
the ubiquity of supermassive black holes in galactic spheroids 
inferred from the M^-a v relation (e.g., Tr emaine et al.ll2002h 
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FIG. 17. — Rest-frame (U-B, Mb) color-magnitude diagram for the objects 
i n the high-S/N sample with [N II]/Ho: emission equivalent width ratio; see 
§ 14.31 The spectral baseline restricts the [N II]/Ha measurements to z < 0.39 
objects. The symbols indicate the equivalent width ratio of Seyfert/LINER 
[log([NlI]/Ha) > -0.3; circle] and Hll-region [log([NlI]/Ha) < -0.3; 
star]; the objects without detectable Ha emission do not have equivalent 
width ratios (square). The colors are as in Fig. [5] and outflows are marked 
by filled symbols. (A color version of this figure is available in the online 
journal.) 



and th eir co-evolution with the lookback time (e.g. JWoo et al.l 
2006) suggests that the energy output from AGNs may play 
important roles in galaxy formation. While AGNs offer in- 
creasingly attractive solutions to the yet elusive mechanism 
for shutting off star formation in massive objects, where and 
how their energy output couples to the gas in and around 
galaxies remains to be identified. Recently, however, sev- 
eral observational studies have elucidated the connection be- 
tween poststarburst and nuclear activity: host galaxies of 
quasars often show poststarburst signatu re in their contin- 
uum emission (e.g.. ICanalizo et ail 120061) : the morphology 
of poststarburst galaxies ofte n show a blue c ore, which 
produces a LINER spectrum dYang et al.l l2006t) : the opti- 
cal emission-line ratios classi fy a majority of poststarburst 
galaxies into Seyfert/LINERs ([ Tremonti et al.l2005l : lYan et al.1 
l200g:lYan & DEEP2 Teamll2006l) . 

Low-ionization outflows are found to be faster in star- 
burst galaxies with some indication of AGN (Ma rtini 120051: 
Rupke et al. 2005b). In z ~ 0.5 poststarburst galaxy sample, 
Tremonti et alj d2007l) find Mg II outflows of ~ 1000 km s~\ 
much faster than the typical Na I D outflows in ULIRGs 
(~ 300-400 km s" 1 ). These fast outflows in the galaxies 
at their postquasar phase provide compelling evidence for 
AGN-driven outflows in poststarbursts. It is worth noting that, 
due to the difference in sample selection, our poststarburst 
outflows in general are l ikely to be of a more typical kind than 
those of iTremonti et al.L i.e., the event which led to the sup- 
pression of star formation does not have to accompany quasar 
activity for us to classify them as poststarbursts. Although the 
exact values are physically meaningless (§ 12. 5b . our measure- 
ments imply Na I D outflow velocities of ~ 100 km s" 1 , which 
do not favor one scenario over others as to the outflow driv- 
ing mechanism. It is plausible that some of our red-sequence 
outflows have g one through the kind of postquasar phase that 
ITremonti et al.l observed, but they do not all have to be. 

Nonetheless, the fact that typical poststarburst galaxies 
show low-level nuclear activity seems convincing, given the 
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FIG. 18. — Rest-frame (U— B, NUV— Rab) color-color diagram for the ob- 
jects in the high-S/N sample wit h [N II]/Ha emission equivalent width ratio. 
Symbol shapes are as in Fig. |17| Symbol colors are as in Fig. [2] and outflows 
are marked by the filled sym bols . The photometric color delineations are de- 
scribed in the caption of Fig. 1131 The spectral baseline restricts the [N II]/Ha 
measurements to z < 0.39 objects. (A color version of this figure is available 
in the online journal.) 

recent studies of the line emission from red-sequence galax- 
ies. In Fig. [17] we show the classification of Seyfert/LINER 
versus star formation using the [Nil] /Ha emission equiva- 
lent width ratio, such that log([N Il]/Ha) = -0.3 divides the 
two classes. 14 The spectral baseline restricts the [Nll]/Ho; 
measurements to z < 0.39 objects. The use of [Olll]/H/3 
allows us to extend the classification to z > 0.39 objects in 
which we see a majority of outflows, but the delineation 
of Seyfert/LINER and star formation becomes notoriously 
ambiguous when o nly that ratio is used. A comparison to 
IWeiner et al.1 (120071) shows that all z > 0.35 outflows in our 
sample have H Il-region-like ratios in [O III] /Hf3, which is 
rather puzzling, given many of them are red galaxies where 
line emissions tend to originate from AGN/LINERs, accord- 
ing to local studies. Using [Nll]/Ho;, Fig. [17] shows ^dis- 
tribution of objects grossly consistent with I Yan et al.1 ( 2006); 
i.e., line emission from the red sequence is dominated by 
Seyfert/LINER, while that from the blue cloud is mostly from 
star formation. Host galaxies of outflows appear to follow the 
parent distribution, yet the analysis suffers from small number 
statistics. 

In Fig. [18] we have a slightly different view on the color 
distribution of galaxies in terms of the emission excitation. 
The advantage of adding the near-UV photometry to detect 
low-level star formation has been discussed in § 13.21 The 
small number of outflow host galaxies at z < 0.39 still lim- 
its our interpretation. Nevertheless, it is quite interesting that 
the transition from star formation to Seyfert/LINER, and then 
to no Ha emission appears to happen along the stellar age se- 
quence, and the abundance of Seyferts/LINERs is observed in 
the region occupied by the galaxies in transition, a significant 
fraction of which appears to host outflows at z ~ 0.5 (Fig.fTSl. 
Although the temporal co-existence of relic outflows and cur- 

14 Unfiuxed DEEP2 spectra force us to use equivalent widths in place of 
line fluxes, yet the adjacent lines forming the ratios are so close that they give 
similar results ( Kobulnicky & Phillips 2003). The Hq-[N II] complex is fitted 
with three Gaussians for emission lines. When Ha absorption is obviously 
present, another Gaussian is fitted; otherwise, a fiducial rest-frame equivalent 
width of 2 ± 1 A is added to the Hct emission flux. 
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rent nuclear activity does not necessarily imply any causal 
connection, extending this analysis to higher redshift, where 
more outflows are expected, could help us understand the po- 
tential impact of winds in establishing the observed trend in 
terms of their driving mechanisms. 

The presence of inflows in red, early-type (Fig. [T3J galax- 
ies in the above figures is also remarkably consistent with 
their host either being quiescent or having their optical emis- 
sion lines excited by AGN/LINERs, with little evidence for 
star formation. It remains to be seen whether any evidence 
exists to directly connect the inflows and the nuclear activi- 
ties. On the other side of outflows, the nature of inflows also 
needs to be explored, since it has been long speculated that 
some mechanism, an AGN being a prime suspect, is prevent- 
ing more stars from being formed than observed in massive 
elliptical galaxies. 

5. SUMMARY 

We reported on a S/N-limited search for low-ionization out- 
flows using the DEEP2 spectra of the 0.11 < z < 0.54 objects 
in the AEGIS survey. Doppler shifts from the host galaxy 
redshifts were systematically searched for in the Na I D op- 
tical resonance absorption doublet. This was the very first 
time that a signature of galactic superwind was systematically 
looked for in the individual galaxy spectra from a modern, 
large spectroscopic redshift survey. 

Our Na I D profile fitting method closely followed that 
of iRupke et alj d2005af) . explicitly fitting the absorption line 
model parameterized by the wavelength and the optical depth 
at the line center, Doppler width, and covering fraction in a 
self-consistent manner. The confidence intervals in the Na I D 
velocities were estimated through the MCMC sampling tech- 
nique. This allowed us to evaluate the quality of Na I D ve- 
locity measurements visually in terms of probability distribu- 
tions of model parameters. Although the spectral resolution 
and S/N limited us to studying the interstellar gas kinemat- 
ics by fitting a single doublet component to each observed 
Na I D profile, LIRG-like outflows should have been detected 
at > 6(7 in absorption equivalent width down to the survey 
limiting S/N (~ 5 pixel" 1 ) in the continuum around Nal D. 
This meant that, if a Na I D o utflow of a compara ble strength 
to the LIRGs detected by the lRupke etal . (2005b) survey left 
its signature in a spectrum, we were able to detect the pres- 
ence in our survey. We discussed the challenges involved in 
recovering physically important parameters from a moderate 
S/N spectrum yet argued that a Doppler shift could be mea- 
sured robustly. 

The detection rate of LIRG-like outflow clearly showed an 
increasing trend with star-forming activity and infrared lumi- 
nosity. However, by virtue of not selecting our sample on star 
formation, we also found a significant fraction of outflows in 
galaxies on the red sequence in the rest-frame (U - B, Mb) 
color-magnitude diagram. Most of these red-sequence out- 
flows were of early-type morphology and showed the sign of 
recent star formation in their UV-optical colors; some showed 
enhanced Balmer H[3 absorption lines indicative of poststar- 
burst as well as high dust extinction. 

We also note that inflows were detected in some red, early- 
type galaxies. Although the definition of an inflow in this 
study was just symmetrically opposite to that of an outflow, 
the fact that we observed them in very distinct subsets of 
galaxies strongly indicated that our outflow/inflow detections 
were not spurious. However, with the difficulty in removing 
stellar absorption component at Na I D as well as the lack of 



immediate explanations for their driving mechanisms, the in- 
vestigation of inflows was beyond the scope of this paper and 
will be explored in future AEGIS studies. We merely note that 
a connection is suspected between an inflow and the feeding 
of nuclear activity in the maintenance mode of AGN feed- 
back. 

The fact that many of our outflows have been found in 
galaxies presumably in transition suggests that galactic su- 
perwinds could outlive starbursts and play a role in quench- 
ing star formation in the host galaxies on their way to the 
red sequence. The detectable imprints of gaseous feedback 
in these galaxies provide us a means to observationally con- 
strain different feedback models. Despite that the small num- 
ber statistics as well as selection effects hindered our ability 
to rigorously characterize the nature of the host galaxies of 
outflows across a wide array of physical parameters accessi- 
ble in the AEGIS survey, the initial analysis presented in this 
paper will help design future experiments on the impact of 
superwinds on galaxy evolution in the epoch when the star- 
forming properties of galaxies drastically change since z <~ 1. 
Gaseous kinematics adds to future studies a useful dimension 
to explore and opens up a promising avenue for constraining 
the driving mechanisms of baryons cycling through the com- 
ponents that constitute the luminous part of the universe, as 
well as for directly quantifying how much gas joins in such 
process. 

T.S. would like to thank the following scientists for inspi- 
rations: Christy Tremonti for sharing her stimulating results 
on SDSS Na I D outflows, as well as her hospitality during 
his visits to Steward; Philip Marshall for enlightenment with 
the Bayesian statistical methods; and David Rupke and his 
coworkers for the series of detailed work on Nal D outflows. 
T.S. gratefully acknowledges Alison Coil for her thorough 
reading of the manuscript and her very insightful suggestions 
and also wishes to extend his thanks to Ben Weiner, Sandy 
Faber, and Francois Schweizer for helpful discussions. We 
thank the anonymous referee for his or her thorough reading, 
constructive feedback, as well as patience. This work was 
never possible without the dedicated efforts, contributions, as 
well as the generosity from all the AEGIS/DEEP2 members. 

The research presented in this paper made an extensive use 
of the Python programming language and the associated tools. 
PyRAF and PyFITS are products of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA. The 
figures in this paper are all prepared by Py Tioga 15 , an open 
source software for creating figures and plots using Python, 
PDF, andTeX. 

This research has made use of the NASA Astrophysics Data 
System abstract service. 

Financial support was provided by the David and Lucille 
Packard Foundation. 

Funding for the DEEP2 survey has been provided by NSF 
grant AST-0071048 and AST-007 1 198. Some of the data pre- 
sented herein were obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the Cal- 
ifornia Institute of Technology, the University of California 
and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial sup- 
port of the W.M. Keck Foundation. 

We gratefully acknowledge NASA's support for construc- 
tion, operation, and science analysis of the GALEX mission, 

15 Py Tioga is available at http://pytioga.sourceforge.net/ 



LOW-IONIZATION OUTFLOWS IN AEGIS 



17 



developed in cooperation with the Centre National d' Etudes Last but not least, we recognize and acknowledge the very 

Spatiales of France and the Korean Ministry of Science and significant cultural role and reverence that the summit of 

Technology. Mauna Kea has always had within the indigenous Hawaiian 

For t he full acknowledg ement of the AEGIS data set, please community. We are most fortunate to have the opportunity to 

refer to lDavis et al.l (120071) . conduct observations from this mountain. 

APPENDIX 

MODELING THE NA I D ABSORPTION LINE 

We closely follow the model presented by Rup ke et ail (|2002, 2005a); readers are strongly encouraged to refer to these papers 
for a thorough discussion. Only a single component of the Na I D doublet is fitted in each spectrum; that is, we assume that our 
fitting of one absorption doublet is sensitive to the bulk propert y of the multiple N a I D absorbing clouds integrated along the 
sightline. While this is certainly an oversimplification (see, e.g.. lRupke et al.ll2.002l for how complex an absorption profile with 
several kinematic components can appear), the limited S/N does not allow more detailed analysis; for the cases in which observed 
Nal D absorption lines are well defined, the fitting results are reliable. For the sole purpose of detecting a Doppler shift from 
a systemic redshift (and not measuring the exact velocity value), a single component fit is an acceptable compromise. We are 
effectively addressing whether an observed Na I D line profile can accommodate a single Doppler-shifted Na I D absorber. 

The observed intensity profile /(A) of an absorption line is fitted by a model profile parameterized in the optical depth space. 
That is, if the (velocity-independent) partial covering fraction is Cf, each Na I D doublet is modeled by 

/(A) = 1 -Cf [1 - e -WA)-T rel] (A)] ^ 

where Tbi ue (A) and T le d(A) are the optical depths of blue and red components of the doublet as a function of wavelength A. (This 
expression is appropriate for a continuum-divided spectrum.) We assume that the velocity distribution of absorbing atoms within 
a cloud is Maxwellian, such that each absorption line is modeled as 

T (A) = T 0e - (A - A ° )2 / (A »W e)2 , 

where to is the optical depth at the line center Ao, c the speed of light, and bo is the Doppler pa rameter in units of speed. Since 
the ratio of oscillator strengths for blue and red sides of the Na I D doublet is two dMortonl 19911) . we may assume that the central 
optical depths are related via ro,biue/2 = ro^ed- Hence the intensity profile of each Na I D doublet component is 



/(A)= 1-C/ jl-exp 



where Abi ue = 5889.9512 A a nd A re d = 58 9 5.9243 A are the rest-frame central wavelengths (in air) of blue and red lines of 
Na I D doublet. As noted by Rup ke etaP (l2005al) . the Maxwellian velocity distribution and velocity-independence of partial 
covering fraction are significant assumptions. In case there is an outflow, an observed Na I D absorption profile likely arises from 
several absorbing clouds entrained in a superwind (e.g., A. Fujita et al. 2008, in preparation), so their bulk kinematics would 
not be described simply by a Maxwellian distribution. Furthermore, if for example a spherical virialized cloud cuts through 
a sightline, the covering fraction must depend on the velocity of the constituent particles in the cloud; thus the assumption of 
velocity independence for covering fraction is not physically consistent in detail. Nevertheless, our analysis does not benefit from 
relaxing these constraints, as the quality of data does not warrant that such detailed information can be extracted. 

Before the above model is fitted, each spectrum around Na I D is divided by the pseudocontinuum, a straight line fitted to 
the variance-weighted spectra within the rest-frame regions of 5822-5842 A and 5910-5930 A; the ranges are chosen after 
visual inspection for the best continuum normalization centered around Na I D, while avoiding other prominent stellar absorption 
features as much as possible. In a very strong outflow, the nebular emission line He I A 5876 can contaminate the bluest wing of a 
Doppler-shifted Na I D component. Visual inspection indicates that few spectra suffer from such a contamination, so no account 
is taken for the He I emission in our measurements. 

We take A re d , bo, To, and Cf as model parameters to be estimated via the Metropolis-Hastings algorithm (Metro polis et alJl953h 
lHastingsl fr970). which yields a more robust confidence interval on a model parameter than that derived from a covariance matrix, 
whe n the probab ility distribution of the model parameter cannot be described as Gaussian. The sampling method also improves 
over lRupke et alj in a sense that each Monte Carlo realization is not drawn from the "best" model which needs to be chosen 
a priori via least-square fitting, for example. The priors for the model parameters are assumed to have a uniform probability 
density over their upper and lower limits: 42 km s -1 < bo < 700 km s -1 , < to < 10 3 , < Cf < 1, and Ao is constrained to be 
within ±700 km s" 1 of the systemic velocity. The \ 2 probability distribution for the given degrees of freedom (i.e., the number 
of data points fitted minus the number of model parameters) is assumed for the likelihood function. Each sampling consists of 
10 5 iterations. 

The measurement pipeline is built on top of PyMC, 16 which implements the Metropolis-Hasting algorithm as an MCMC 
sampler, and developed in Python. The distributions of model parameters are visually inspected with a graphical user interface 
(GUI) along various dimensions; see Fig. [19] The integrity of fitting to the data spectrum is checked at several interactively picked 
points on the two-dimensional intensity map of the distributions of model parameters. The GUI-driven visual inspection guards 
against the fits latching on to low-S/N features and helps to identify unphysical fit results. The distribution of A le d, marginalized 
over all the other parameters and convolved with the redshift uncert ainty , needs to be roughly Gaussian and well bounded within 
±700 km s" 1 to make it into the high-S/N Na I D velocity sample (§ 12.41 ). The software may be open sourced at a later date. 
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FIG. 19. — Few output windows from the Nal D measurement pipeline. In this figure, the distributions of model parameters are shown as the two-dimensional 
intensity color maps. Any location on the graphical user interface (left) can be clicked on to extract the model spectrum overlayed on top of the data spectrum 
(right). The relatively high-(S/N) p j x spectrum yields well-behaving parameter distributions here. The figure also gives an example of typical model parameter 
distributions. On the left, the parameters plotted are To-A re d (top left), foo-A rel i (top right), C/--A re( i (bottom left), and Cf-rg (bottom right). If marginalized over 
the parameter against which they are plotted, the distributions of A lc( j and bo become roughly Gaussian, where as those of Cf and to are not. It is apparent that 
Cf is not well constrained for Cf > 0.45, and to, while relatively well constrained, is highly correlated with Cf. See § |2.4l for detail. (A color version of this 
figure is available in the online journal.) 
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